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Abstract 

A brief introduction on the issue of stability in generalized modified gravity is presented and the 
dynamical system methods are used in the investigation of the stability of spatially flat homogeneous 
cosmologies within a large class of generalized modified gravity models in the presence of a relativistic 
matter-radiation fluid. 
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O ■ 1 Introduction 

<N 

■ To start with, we recall that recent cosmological data support the fact that there is a good evidence for 

a late accelerated expansion of the observable universe, apparently due to the presence of an effective 
j /^ ■ positive and small cosmological constant of unknown origin. This is known as dark energy issue (see for 

example [I]). 

Among several existing explanation, the so called modified gravity models are possible realizations 
of dark energy (for a recent review and alternative approaches see [2j [3j |4]), which may offer a quite 
natural geometrical approach again in the spirit of the original Einstein theory of gravitation. In fact, 
the main idea underlying these approaches to dark energy puzzle is quite simple and consists in adding to 
the gravitational Einstein-Hilbert action other gravitational terms which may dominate the cosmological 
evolution during the very early or the very late universe epochs, but in such a way that General Relativity 
remains valid at intermediate epochs and also at non cosmological scales. 

The A-CDM model is the simplest possibility but, it is worth investigating more general modifications, 
possible motivations run from quantum corrections to string models. We shall first consider the simpler 
modification of the kind F(R) = R + f{R), and then discuss the related generalizations. Models based 
on F(R) are not new and they have been used in the past by many authors, for example as models for 
inflation, f(R) = aR 2 [5 . Recently their interest in cosmology was triggered by the model f(R) = —fi 4 /R, 
proposed in order to describe the current acceleration of the observable universe [51 [7] . For incomplete 
list of references, see [5]. 

It is important to stress that these F(R) models are conformally equivalent to Einstein's gravity 
coupled with a self-interacting scalar field, Einstein frame formulation. We will consider only the Jordan 
frame, in which the dynamics of gravity is described by F(R) with minimally coupled matter. Observa- 
tions are typically interpreted in this Jordan frame. 

We would like also to mention the so called viable F(R) models, which have recently been proposed 
Uni HU O EH [H] , with the aim to describe the current acceleration with a suitable choice of F(R) — 
R + f(R), but also to be compatible with local stringent gravitational tests of Einstein gravity F(R) = R. 
The main idea is the so called disappearing of cosmological constant for low curvature, and mimicking 
the A-CDM model for high curvature. Thus, the requirements are: 

a. f(R) — * , R — ► , compatibility with local tests; 

b. f(R) — ► — 2Ao , R — ► +oo, description of current acceleration; 

c. Local stability of the matter. 

As a illustration, we recall a recent example of viable model [IS] 

f(R) = -a (tanh ^ R ~^ + tanh ^ 

where Rq, and Ao are suitable constants. Its advantages are a better formulation in the Einstein frame 
and a generalization that may also include the inflation era. 
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2 The de Sitter stability issue for F(R) models 



The stability of the de Sitter solution, relevant for Dark energy, may be investigated in the F(R) models 
in several ways. We limit ourselves to the following three approaches: 

i. perturbation of the equations of motion in the Jordan frame; 

ii. one-loop gravity calculation around de Sitter background; 
hi. dynamical system approach in FRW space-time. 

We shall briefly discuss first two approaches, and then concentrate to the third one, since it is the 
only that can be easily extended to more general modified gravitational models, which is the main issue 
of this short review. 

i. Stability of F(R) model in the Jordan frame 

The starting point is the trace of the equations of motion, which is trivial in Einstein gravity R = —k 2 T, 
but, for a general F(R) model, reads 

3V 2 f'(R) - 2f(R) + Rf'(R) — R = n 2 T . 

The new non trivial extra degree of freedom is the scalaron: 1 + f'(R) — e~ x . Requiring R = Rq = CST, 
one has de Sitter existence condition in vacuum 

R o + 2f(Ro)-R o f'(R o ) = 0. 

Perturbing around dS: R = Rq + SR, with SR = — ^prj|~) &Xi one arrives at the scalaron perturbation 
equation 



V 2 8 X ~ M 2 S X 



K 
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6(1 + f'(R )Y 
One may read off the scalaron effective mass 



3 V f"(Ro) 

Thus, if M 2 > 0, one has stability of the dS solution and the related condition reads 



Rof"(Ro) 



> 1. 



If M 2 < 0, there is a tachyon and instability. Furthermore, one may show that M 2 has to be very 
large in order to pass both the local and the astronomical tests and 1 + f'(R) > 0, in order to have a 
positive effective Newton constant. The same result has been obtained within a different more general 
perturbation approach in [16) . 



ii. One-loop F(R) quantum gravity partition function 

Here we present the generalization to the modified gravitational case of the study of Fradkin and Tseytlin 
[17j . concerning Einstein gravity on dS space. One works in the Euclidean path integral formulation, 
with dS existence condition 2^0 = RqFq, assumed to be satisfied. The small fluctuations around this dS 
instanton may be written as 

9ij = 9o,ij + hi , .'/'•' = 9% - /''' + h ik h{ + O (h 3 ) , h = 9% hj . 

Making use of the standard expansion of the tensor field hij in irreducible components, and making an 
expansion up to second order in all the fields, one arrives at a very complicated Lagrangian density L2, 
not reported here, describing Gaussian fluctuations around dS space. As usual, in order to quantise 
the model described by L2, one has to add gauge fixing and ghost contributions. Then, the computa- 
tion of Euclidean one-loop partition function reduces to the computations of functional determinants. 
These functional determinants are divergent and may be regularized by the well known zeta-function 
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regularization. The evaluation requires a complicated calculation and, neglecting the so called multi- 
plicative anomaly, potentially present in zeta- function regularized determinants (see j!8j). one arrives at 
the one-loop effective action [H] (here written in the Landau gauge) 



247rF f 



on — shell — 
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The last term is absent in the Einstein theory. As a result, in the scalar sector one has an effective mass 
— Rq ) . Stability requires M 2 > 0, in agreement with the previous scalaron analysis, and 



M 2 = 



1 / 2F 

3 \WF^ 



with the inhomogeneous perturbation analysis [16j . 



The dynamical system approach 

The main idea of such a method is to convert the generalized Einstein- Friedman equations in an au- 
tonomous system of first order differential equations and makes use of the theory of dynamical systems 
(see [201 [2"T1 [2"2l [231 [2~4l [2"5] and references therein) . We remind that the stability or instability issue is 
really relevant in cosmology. For example, in the ACDM model it ensures that no future singularities 
will be present in the solution. Within cosmological models, the stability or instability around a solution 
is of interest at early and also at late times. 

The stability of de Sitter solutions has been discussed in several places, an incomplete list being 
[H1[M1[571 ?, [IE [IS [101 CHI Gill 111]- More complicated is the analysis of other critical points, associated 
with the presence of non vanishing matter and radiation. For example, in the paper [34] , a non local 
model of modified gravity F(R) has been investigated by means of this approach. 

Here we shall extend to general modified gravity in the presence of matter the method given in 
Ref. [22] . which permits to determine all critical points of a F(R) model. Ordinary matter is important 
in reconstructing the expansion history of the Universe and probing the phenomenological relevance of 
the models (see for example the recent papers |22[ [2"3l [2"4l [2"5] . where the F(R) case has been discussed 
in detail). Our generalisation consists in the extension of that method in order to include all possible 
geometrical invariants. This means that F could be a generic scalar function of curvature, Ricci and 
Riemann tensors. 

There are some theoretical (quantum effects and string-inspired) motivations in order to investigate 
gravitational models depending on higher-order invariants. The "string-inspired" scalar-Gauss-Bonnet 
gravity case F(R, G) has been suggested in Ref. pS] as a model for gravitational dark energy, while 
some time ago it has been proposed as a possible solution of the initial singularity problem 36J. The 
investigation of different regimes of cosmic acceleration in such gravity models has been carried out in 
Refs. [HZl [Ml EH1 HOI IHl IHl H31 IIH HH1 SS]. In particular, in 03] a first attempt to the study of the 
stability of such kind of models has been carried out using an approach based on quantum field theory. 

The method we shall use in the present paper is based on a classical Lagrangian formalism, see, for 
example, (47[ I48[ I49j , inspired by the paper [5] , where quantum gravitational effects were considered 
for the first time. With regard to this, it is well known that one-loop and two-loops quantum effects 
induce higher derivative gravitational terms in the effective gravitational Lagrangian. Instability due to 
quadratic terms have been investigated in [50j . A particular case has been recently studied in [51j and 
general models depending on quadratic invariants have been investigated in [521 153j . 

A stability analysis of nontrivial vacua in a general class of higher-derivative theories of gravitation 
has already been presented in |54j . Our approach is different from the one presented there since we are 
dealing with scalar quantities and moreover it is more general, since it is not restricted to the vacuum 
invariant submanifold. 

Finally, it should be stressed that the stability studied here is the one with respect to homogeneous 
perturbations. For the F(R) case, the stability criterion for homogeneous perturbations is equivalent to 
the inhomogeneous one [16j . In the following we will summarize the results obtained in |55j . 
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2.1 The dynamical system approach: The general case 

To start with, let us consider a Lagrangian density which is an arbitrary function of all algebraic invariants 
built up with the Riemann tensor of the FRW space-time we are dealing with, that is 

£=-^-F(R,P,Q,...)+C m , (2.1) 

where R is the scalar curvature, P = R^R^ and Q = R al3l5 R a p-yS are the two quadratic invariants and 
the dots means other independent algebraic invariants of higher order, and is the matter Lagrangian 
which depends on p and p = p(rho), the density and pressure of the matter. 
For the sake of convenience we write the metric in the form 

ds 2 = -e 2n{t) dt 2 + e 2a{t] dx 2 , N(t) = e" (t) , a{t) = e a(t) . (2.2) 

In this way a(t) — H(t) is the Hubble parameter and a generic invariant geometrical quantity U has the 
form 

U = e~ 2pn ^ u(h, a, a) = e ~ 2pn ^ u(h, H, H) = H 2p e " 2p " (t) u(X) , (2.3) 

where 2p is the dimension (in mass) of the invariant under consideration and X = (H / H 2 — h/H) (see 
55 1). In particular one has 

R = 6H 2 e- 2n (2 + X), 

P=12H i e~ 4n {3 + 3X + X 2 ), (2.4) 
Q = 12H 4 e- 4n (2 + 2X + X 2 ) , 

Using this notation, the action reads 

S = - J d 3 x J dtL(n,n,a,a,6t) = ^- J d 3 x J dte n+3a F(n, h, a, a) + S m , (2.5) 
and the Lagrange equations corresponding to the two Lagrangian variables n(t) and a(t) are given by 

En = ^ - |§ = 2^T 005 00 = 2p , (2.6) 

dL d dL d 2 dL . m nh n . 

oa dt aa dt z act 

It has to be noted that since n(t) is a "gauge function" , which corresponds to the choice of the evolution 
parameter, Eqs. (|2.6p and (|2.7p are not independent and in fact they satisfy the differential equation 

^ = hE n + aE n , (2.8) 

which is equivalent to the conservation of energy-momentum tensor. Furthermore, we may use the gauge 
freedom and fix the cosmological time by means of the condition N(t) = 1, that is n(t) = 0. From now on 
it is understood that all quantities will be evaluated in such a gauge and so the parameter t corresponds 
to the standard cosmological time. In this way Eqs. (|2.6|) and (|2.7|) read 

HF k - HF H + F- HFfj + 3i? 2 F^ = 2p , (2.9) 

F ti -F H + SHFfj - 3HF H + 3F + 3HF^ + 9H 2 F^ = -6p . (2.10) 

The latter equations are the generalisation to arbitrary action of the well known Friedmann equations. 

Now we shall replace Eqs. (|2.9p - (|2.10p with an autonomous system of first order differential equations. 
To this aim we first observe that in pure Einstein gravity, that is for F = R, (f2~9l) - (|2~T0|) read 

H 2 F tl = F x = 2p => fl p = l, (2.11) 

(3H 2 + 2H)F ti = (3 + 2X)F X = -Gp flp = -l-~X , (2.12) 

o 
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where we have introduced the dimensionless variables 

ttp = ^ = ^, tt p = ^ = ^, (2.13) 

9 E 2 F ti F x P H 2 F ti F x 

which in this special case are given by the usual values f2 p = p/3H 2 and £l p = p/3H 2 . From Eqs. (|2.11| 
and (j2~T2")l it follows 

(2.14) 

p Hp 6 

In the general case, Eqs. (|2.9|) and (|2 . 10[> have more terms with respect to (|2.11| and (12.121) and it is 
quite natural to interpret them as corrections due to the presence of higher-order terms in the action. 
Then we define 

n e /f = n p + n c p urv = i , f^ = n p + n£«™ = -i-|x, (2.15) 

Qeff 2 

where fi™ru and Sip"™ are complicated expressions, which only depend on the function F. They can be 
derived from (|2.9|) and (|2.10|) . but they explicit form is not necessary for our aims. The effective quantity 
w e f f is equal to the ratio between the effective density and the effective pressure and it could be negative 
even if one considers only ordinary matter. It is known that the current-measured value of w e ff is near 
-1. 

In order to get all critical points of the system now we follow the method, as described, for example, 
in [22] • First of all, we introduce the dimensionless variables 

Q p = ^ = ^, <7 P =J^=^, (2.17) 

9 E 2 F ti F x P H 2 F k F x 

X=*L, Y = F - 9 HFh = — -X, Z= Fh ~ Fh =Ek-2X-Z, (2.18) 
H 2 ' H 2 F ti F x HF h F x 

where the prime means derivative with respect to a and the quantity 

has to be considered as a function of the variables X and Y. In general it is a function of X and H , but 
this latter quantity can be expressed in terms of X and Y as a direct consequence of the definition of Y 
itself. Then we derive an autonomous system by taking the derivatives of such variables. From Eq. (|2.9[) 
we have the constraint 

n p = Y + Z + 3 =!> il curv = -(Y + Z + 2), (2.20) 

which reduces to the standard one when F is linear in R (general relativity with cosmological constant). 

Deriving the variables above by taking into account of (|2.10p and (I2.20|l we get the system of first 
order differential equations 

X' = -2X 2 - jX + /3(Z + f) 

Y' = -{2X + Z + i)Y-XZ (2.21) 
Z' = -3(1 + w){Z + Y + 3) - (Z + £){Z + 3) - X(Z + 6) 

where X' = ^2L = (and so on) and for simplicity we have set p = wp, with constant w. For ordinary 

matter 0<w<l/5(w — corresponds to dust, while w = 1/3 to pure radiation), but in principle one 
could also consider "exotic" matter with w < or cosmological constant which corresponds to w = — 1. 
We have also set 

= P(X,Y) = - F S— = -pL, 7 = 7i x,Y) = - F f?- = ?p^=ft x +Z. (2.22) 
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It is understood that Fjrjj =/= has been assumed. The quantity f2 p at the critical points will be 
determined by means of Eq. (|2.20[) . 

The critical points are obtained by putting X' = 0,Y' = 0, Z' = in the system (|2.2ip . The number 
and the position of such points depends on the Lagrangian throughout the functions /3, 7 and £. In 
principle, given F one can derive all critical points, but in practice for a generic F the algebraic system 
could be very complicated and the solutions quite involved. We shall consider in detail some particular 
cases at the end of the Section. 

As already said above, the critical points of (|2.2ip depends on £, (3, 7, which in general are complicated 
functions of X and Y, then it is not possible to determine general solutions without to choose the model, 
nevertheless it is convenient to distinguish two distinct classes of solutions characterised by the values of 
w =/= —1 and w = — 1. For the sake of completeness we consider w < 1/3 and so we write the solutions 
also for "exotic" matter, that is quintessence (—1 < w < 0) and phantom (w < —1). Of course, such 
solutions have to be dropped if one is only interested in ordinary matter/radiation. We have 

• w — 1 — The critical points are the solutions of the system of three equations 
2X 2 + jX - (3(Z + = 

(2X + Z + £)Y + XZ = (2.23) 
3(1 + w)(Z + Y + 3) + (Z + £){Z + 3) + X(Z + 6) = 

where £, f3, 7 are functions of X, Y determined by Eqs. (|2. 19|) and (|2.22j) . The stability matrix has 
three eigenvalues and the point is stable if the real parts of all of them are negative. 

The latter system has always the de Sitter solution Pq = (X = 0, Y = 1, Z = —4), where fl p = 
and w e ff = — 1. Note however that such a solution could exist also in the presence of matter, since 
the existence of Pq critical point only implies that the critical value for Q p vanishes. 

f2 p =/= 0, w = — 1 — The critical points are given by 
2X 2 +jX-(3(Z + = 

(2X + Z + £)Y + XZ = w = -l. (2.24) 

(Z + 0(Z + 3) +X(Z + 6) = 

For this class of solutions, the non-singular stability matrix has three eigenvalues and the point is 
stable if the real parts of all of them are negative. 

We see that there is at least one singular case (critical line) when X = and Z = — £ = —4. In fact 
in such a case Y or Q p are undetermined since 

n p =Y + 3- £(0, Y) = Y - 1 => Y = 1 + n p , fl p arbitrary . (2.25) 

Such a solution can be seen as a generalisation of the de Sitter solution for a model with cosmological 
constant. The de Sitter critical point for the model F = F — 2A reads (X = 0, Y = 1,Z = —4). 
Such a solution follows from Eq. (12.25)) if we choose po = A. In fact, on the critical point (X = 
0, Y = 1 + n p , Z = -4) (Eq. (|2~2"51) ) and from definitions (|2~Hj|) we get 

"^wk'WF-- 1 * fs ^k = U (2 - 26) 

which corresponds to de Sitter critical point for F. Eq. (12.25j) is more general than the case with 
pure cosmological constant since p is not necessary a constant, and for this special class of solutions 
w e f f = — 1. Note also that the stability matrix has always a vanishing eigenvalue and the stability 
of the system is determined by the other two eigenvalues. 

For some models, but just for technical reasons, it could be convenient to treat the cosmological 
constant as matter, using the previous identification we have done. 
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Explicit examples 

In order to see how the method works, now we give explicit solutions for some models and, when possible 
we also study the stability of the critical points. We restrict our analysis to the values < w < 1/3 and 
to the special value w — —1, which corresponds to the pure cosmological constant, but in principle any 
negative value of w could be considered, even if this will be in contrast with the aim of modified gravity. 
In fact, modified gravity can generate an effective negative value of w without the use of phantom or 
quintessence. 

It as to be stressed that in general, due to technical difficulties, one has to study the models by a 
numerical analysis. Only for some special cases one is able to find analytical results. Here we report the 
results for some models of the latter class in which the analytical analysis can be completely carried out. 
We also study more complicated models and for those we limit our analysis to the de Sitter solutions. 

In the following we shall use the compact notation 

P=(X,Y,Z,Q p ,w eff ), P = (0, 1,-4, 0,-1), P A = (0,l + fi A ,-4,ftA,-l). (2.27) 

The latter is an additional critical point that we have for the choice w = — 1 and can be seen as the de 
Sitter solution in the presence of cosmological constant. 

F = R — /i 4 / 1 R — This is the well known model introduced in [7] and discussed in [22] . For this 
model the system (|2.21j) with arbitrary w has six different solutions, but only two of them effectively 
correspond to physical critical points, if < w < 1/3. In principle there are other critical points 
for negative values of w (phantom or quintessence) and moreover there is also a particular solution 
for w = — 1 which corresponds to the model with a cosmological constant A. 

Solving the autonomous system one finds 

• P = Pq: unstable de Sitter critical point. The critical value for the scalar curvature reads 
Ro = V3> 2 - 

• P = (—1/2,-1,-2,0,-2/3): stable critical point. At the critical value, Ho = 0. 

• P = (3(1 + w)/2, -(5 + 3u>), -2(5 + 3u>), -3(4 + 3w), -(2 + w)): unstable critical point where 
H = 0. 

• P = P\: unstable critical point. At the critical value one has Hq — (A/6)(l + y/l-h 3/i 4 /4A 2 . 

F = R + aR 2 + bP + cQ — (Starobinsky-like model). Here we have to assume 3a + b + c ^ otherwise 
the quadratic term becomes proportional to the Gauss-Bonnet invariant. For < w < 1/3, this 
model has only one critical point. In order to have a de Sitter solution, we have to introduce a 
cosmological constant A. We have in fact 

• P = Pq: Minkowskian solution with Rq — 0, which is stable if 3a + b + c > 0. 

• P = P\: de Sitter critical point with Ro — 6A, which is stable if 3a + b + c > 0, in agreement 
with [53] . 

F = R — cPQs — This is the simplest toy model with the cubic invariant Qz — R a ^ s R a @ij,i>Rttg ■ For 
this model we have the following critical points: 

• P = Pq: unstable de Sitter solution with Ro = 6/d. 

• P = Pq: stable Minkowskian solution with Ro = 0. 

• P = (0.05, 0.60, -3.60, 0,-1.03): stable solution with H = 0. 

• P = P\: this point exists, but not for any value of d and A. Also the value of Ho and the 
stability depend on the parameters. 

F = R + aR 2 + bP + cQ — cPQs — This is a generalisation of the previous two models. It may be 
motivated by the two- loop corrections in quantum gravity [561 157] . The de Sitter critical points 
have been studied in Ref. 58 J. The algebraic equations (|2.21|) arc too complicated to be solved 
analytically, but it is easy to verify that there are at least the following solutions: 

• P = Po'. de Sitter solution with Ro = 6/d. This is stable if 3a + b + c + 3d > 0. 

• P = Po: Minkowskian solution with Ro = 0, which is stable if 3a + b + c > 0. 

• P = P\: also in this case this point exists and is stable depending on the parameters (see 

[SHI). 
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3 Conclusion 



In this contribution, after a short review on the F(R) case, we have presented the dynamical system 
approach which has permitted to arrive at a first order autonomous system of differential equations 
classically equivalent to the equations of motion for models of generalized modified gravity based on an 
arbitrary function F(R, P, Q, Q3...), namely built up with all possible geometric invariant quantities of 
the FRW space-time. We have shown that, in the special case of F(R) theories, the method gives rise to 
the well known results [551 H31 HI] > but in principle it can be applied to the study of much more general 
cases. 

As illustrations, we have discussed some simple models, for which a complete analytical analysis 
concerning the critical points has been carried out. However, in general, due to technical difficulties, a 
numerical analysis is required. Among the models investigated, we would like to remind that we were 
able to deal with one which involves a cubic invariant in the curvature tensor and, to our knowledge, this 
has never been considered before, and this shows the power of the present approach. 
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